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Abstract 



^ , The Cross-Newell equations for hexagons and triangles are de- 



rived for general real gradient systems, and are found to be in flux- 
divergence form. Specific examples of complex governing equations 



/\ ' that give rise to hexagons and triangles and which have Lyapunov 



functionals are also considered, and explicit forms of the Cross-Newell 
equations are found in these cases. The general nongradient case is 
also discussed; in contrast with the gradient case, the equations are 
not flux-divergent. In all cases, the phase stability boundaries and 
modes of instability for general distorted hexagons and triangles can 
be recovered from the Cross-Newell equations. 



1 Introduction 

Hexagons are a very common planform arising in pattern-forming systems. 
The asymmetry between the centres and the edges of the hexagons leads 
to a favouring of hexagonal patterns in situations where there is intrinsic 
asymmetry, such as in Benard-Marangoni convection IQ where the top 
surface of the convecting layer is free and the bottom surface is in contact with 
a rigid boundary. Most natural systems will have some degree of asymmetry, 
and hence hexagons are widely observed, not only in convection experiments, 
but also for example in vibrated granular layers 0] and during directional 
solidification 0. Triangular patterns are more unusual, but are seen in some 
systems, such as vibrated granular layers |^. 

Cross and Newell pioneered a method of describing the behaviour of a 
fully nonlinear roll pattern in an extended system by following the evolution 
of the local phase, and hence the wavevector, associated with the roll pattern 
as it varies in space and time. This method was further developed by Passot 
and Newell who regularised the Cross- Newell equations outside the region 
of roll stability, introducing an order parameter equation to account correctly 
for the behaviour of the pattern in regions where the amplitude is small. 

The purpose of the current paper is to apply ideas similar to those of 
Cross & Newell [0 and Passot & Newell [^ to the evolution of fully nonlinear 



hexagonal and triangular patterns in large aspect ratio systems, such as 
those seen in experiments on Rayleigh-Benard convection in SFg near the 
thermodynamical critical point 0. 

The paper is structured as follows: §2 presents a method of deriving the 
Cross-Newell equations for triangles and hexagons in a general real gradient 
system. The Cross-Newell equations for particular complex gradient systems 
are derived in §3 for hexagons and §4 for triangles. The case of free hexagons 
and triangles is discussed in §5, and the general nongradient case in §6. 
Section 7 concludes and indicates some directions for future investigation. 

2 Derivation of the Cross-Newell equations 

It is assumed that fully developed hexagons or triangles can be described by 
a stationary solution w = Wo(x) of an equation Wt = Cw + TVw, where C 
and M are linear and nonlinear operators respectively, at least one of which 
is differential, with variational structure such that 

/— /-^^. (^) 

where w and G{w) are real. 

Hexagons and triangles are described by three wavevectors ki , k2 and ka 
forming a resonant triad such that ki + k2 + ks = 0. In the case where the 



governing equations force this resonance to be maintained, the pattern can be 
described using two phases 6i and 62 associated with two of these wavevectors 
ki = V^i and k2 = V^2- For fully nonlinear triangles and hexagons, the 
hexagon amplitude and the total hexagon phase a = 61 + 62 + 03 (where 
ks = — ki — k2 = V^s) are determined adiabatically from the two phases 61 
and 62 except in the vicinity of defects where the amplitude is small or when 
the driving stress parameter of the system is close to the critical value for 
pattern formation so that the amplitude is small everywhere. In the 'free' 
case discussed in §5, where the resonant triad may be broken, the hexagon 
amplitude is slaved to the three independent phases 61, 62 and 63, except 
when the amplitude is small. 

In a large aspect ratio system, the size and orientation of the hexagons 
will typically change slowly in space and time. To describe these changes, it 
is convenient to introduce large scale phases 9i = e9i, ©2 = e^2, where e -C 1 
is the inverse aspect ratio of the box, and slow space and time scales X = ex, 
T = e^t. The local wavevectors are then given by ki = Vx^i = Vx0i, 
i = 1,2. 

The hexagon solution is now considered to be a function of the two phases 
61 and 62 and the slow space and time scales, such that w = 10(61,62;^, T). 



Hence the space and time derivatives of w are given by 

Vxw(ei,^2;X,T) = (kia,i + k29,2 + eVx)w^(^i,^2;X,T), (2) 

9tw(^i,^2;X,T) = {eQiTdei + eQ2Tde2 + e''dT)w{e,,e2;X,T). (3) 

To leading order then the following equations hold 

Wt = e{eiT(<9eiWo) + 02t(<992Wo)}, (4) 

5w = {d0iwo)59i + {d02Wo)592. (5) 

Substituting all this information into the governing equation ( |^) and aver- 
aging over 6i and 6*2 gives, to leading order in e, 

e j dxdy{QiT{dgiWo) + <d2T{de2Wo)}{ideiWo)60i + {802100)662} 
f dG dG dG 



where (.) denotes the averaging. Remarking that 5k\ = 5(Vx^i-Vx^i) = 
2V^6i.V^66i = 2ki.V^66i, and similarly that (5A;^ = 2k2.Vx(56'2 and(5(ki.k2) = 
ki.Vx(56'2 + k2.Vx(56'i, the divergence theorem can be used with suitable 
boundary conditions, to show that 

e J dxdy{QiT{deiWo) + 02T{de2Wo)}{{deiWo)Sei + {802100)662} 



+eJdxdyV. |(ki5^2 + k2'5^i)^7^^| , (7) 

where V = Vx = e^"'^(Vx — ki^ei — 1^2992) ■ Since 69i and 662 are arbitrary, 
it is possible to extract the phase equations 

dG . dG \ .„. 



Qirldeiwol^ + Q2T{deiWo){d02Wo) = V. 2ki— ^ + k2 



dki 9(ki.k2) 



02T\de2Wo\' + QiT{deiWo){d02Wo) = V. ( 2k27^ + ki ^^^^ ^ J . (9) 



dG dG 



The phase stabihty boundaries for general distorted hexagons and trian- 
gles defined by wavevectors ki and k2 can be recovered from the Cross-Newell 
equations by first writing the equations explicitly in terms of the phases to 
give 



©iTl^ei^op + Q>2T{.deiWo){de2Wo) = 

V. |2vei^(|vei|2, |ve2p, (vei.ves))! + 



V. {v02^^^^(|veip, ive^p, (vei.ve^))! , (lo) 



&2T\de2Wo\^ + &iT{deiWo){d02Wo) = 

V. |2V02^(|vei|2, |ve2p, (vei.ves))! + 



V. {^0i^(^^(iv0il'' iv02p, (vei.ve^))! , (ii) 



and then setting Oi = ki.X + Oi and 62 = k2.X + 02, where 9i and 
B2 are smalL Linearising in 0i and B2, and setting 0i = @^e'^'r+«kx ^^^ 

©2 = ^2^'^^'^^^''^ , with Gi and 62 real constants gives a dispersion relation 
for the growth-rate eigenvalues a. Hence the stability boundaries and modes 



of instability can be found as in ||T0| [|TT|]. Direct numerical integration of the 



Cross-Newell equations could also be used to determine the region of stable 
hexagons and triangles, and comparison could be made with the stability 



region for regular hexagons found by other numerical methods as in p2 . 

Ideally the governing equations should be real, as assumed here, in order 
to allow the formation of disclinations on an individual set of rolls |^. How- 
ever, there do not appear to be simple examples of real governing equations 
which give fully nonlinear hexagons or triangles as an exact stationary solu- 
tion, and so in order to make further explicit analytical progress we shift our 
attention in the following section to complex governing equations which do 
indeed give hexagons. This is perhaps less of a handicap than it would be in 
the case of rolls, since the canonical hepta-penta defect of hexagons is made 
up of dislocations, which can be described by a complex order parameter. 



3 Cross-Newell equations for hexagons 

With slight modifications to the spatial derivative terms, the standard com- 



plex amphtude equations for hexagons ||T3[ [|14| can be used as the basic 
governing equations, giving 



-^ = ^Zi + az*j^^z*j^2 - f^\zi\^Zi - 7(|zi+i|^ + \zi+2\'^)zi + V'^Zi, (12) 



where A, a, /5 and 7 are real constants, and where * denotes complex con- 
jugation. The hexagon solutions are represented hy w = Re{zi + Z2 + Z3), 
Zi = Rie'^^% with i = 1,2,3 and cyclic. Here the usual spatial derivatives have 
been replaced by V^ in order to preserve the isotropy of the system. 

There is a Lyapunov functional associated with the amplitude equations, 
given by 

L = -A(|zi|2 + |z2|' + ksH - a{ziZ2Z3 + z^z^zl) + -/5(|zi|^ + Wf + l^sl^) 

+7(kirk2r + \z2Wz? + \zzWxf) + ivzip + \vz2f + iv^sr, (is) 

such that 

^ = -i^. (14) 

dt 5zt ^ ' 

There are wavevectors associated with the phases according to ki = V6'j as 
before. A hexagonal or triangular pattern arises when the sum of the three 
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wavevectors is zero, i.e. ki + k2 + ks = 0. Hence the total phase X]j ^j = (^{t) 
is a function of time only. 

The fully nonlinear hexagonal solution takes the form w = Ri cos 6i + 
i?2 cos 02 + -R3 cos ^3 where 

= i?i(A - kl) + aRiRs cos a - pRf - -f{Rl + Rl)Ru (15) 

= i?2(A - kl) + aR^Ri cos a - PRl - -f{Rl + Rl)R2, (16) 

= i?3(A - kl) + aRiR2 cos a - (3Rl - -f{Rl + Rl)Rs, (17) 

= asina (18) 

hold, and where the Ri are nonzero constants. Clearly if a is nonzero, as 
assumed in this section, the total phase a must take the value or tt. If a is 
zero, a can take any value. 

In the case of nonzero a, there are only two independent phases, which 
without loss of generality are taken to be 61 and 62- The third phase 6^ = 
a — 9i ^ 62 is then dependent, since a is fixed. 

As in the previous section, it is assumed that the wavevectors vary slowly 
in space and time, so that it is possible to define large scale phases Oj = e6i 
and long space and time scales such that Vx = ^idg^ +k2(902 +eVx and dt = 
e^Qirdg^ + Q2Td02)- The solution is expanded in the form z = (^1,2:2,-23) = 
Zq + ezi + e^Z2 + ..., where Zq is the fully nonlinear hexagon solution above. 

9 



To leading order, the Lyapunov functional for the fully nonlinear hexagons 
takes the form 

L = -X{Rl + Rl + Rl) - 2aRiR2R3 cosa + -{R^ + Ri + RJ) 

+^{R\Rl + RlRl + RlRl) + klRl + klRl + klR^ (19) 

Since ki + k2 + ks = holds, A;| can be rewritten kl + k"^ + 2ki.k2, and it 
is clear that the i?j, a and L all depend only on fc^, k^ and ki.k2. Hence the 
variation 5L in the Lyapunov functional is given by 






It is also clear that 



'' - i^< + 'i'- - ->■ - '4'- (^^ 



holds. Further, it can be seen that 



5z* = {de,z*o)69, + {de,z*o)S92, (22) 

dtz = eOMde^zo) + ee2Tide,zo) + ©(e^) (23) 

hold. Substituting these into equation ( |21|) gives 

5L = -e{Q,T{de,zo) + e2T{de,zo)}.{{dg,z*^)60, + {de^z^SO^} + c.c. (24) 
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to leading order, where c.c. denotes complex conjugate. Considering / SLdxdy, 
where the integral is taken over the whole domain, gives 



dxdye{eiT{dg^zo) + Q2T{de2Zo)}.{idg,zl)Sei + {89^21)662} + c.c. 
f dL dL dL 



2, 



(25) 
from which it is easy to identify the phase equations 



01T(5eiZo).(9,,Z*) + 02T(5e,Zo).(9,,Z*) + C.C. = V. 2ki-^ + k2 



dkl '^(ki.ka);' 
(26) 

e2T(5,,Zo).(9,,Z*)+0iT(5,,Zo).(9,,Z*)+C.C. = V. (2k2^ + ki— — — j. 

(27) 
The straightforward substitutions 



/ i?ie^^i 



zo 



R2e''' I , (28) 
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dejZo = 



(29) 



992^0 



\ 



(30) 



reduce the phase equations to 



2{Rl + Rl)QiT + 2RIQ2T = V. 2ki 



aL 



aL 



aA;2 "^(ki.ka) 



(31) 



2(i.^ + i.De.. + 2i.^e..^V.(2k.g + k.^ 



(32) 



Substituting equations ( |T5])-( |T8|) into equation ( |T9[) gives a simpler expres- 
sion for the Lyapunov functional 



L = aRiR2R3 cos a- -i{R\Rl + RIrI + RIrI) - -I3{R\ + R2 + Rt), (33) 



which when differentiated with respect to k\ becomes 



dL 



'dRi 



dRo 



dR. 



a cos a \—^R2R3 + Ri-Kr2R3 + RiR2-Kp 



dRi 



dR 



2/d2 



dR' 



3/d2 



-27 Ri^m + Ri) + R2^{Ri + Ri) + R^i^m + ^^ 



9A;i^ 



9A;. 



dB 



-2/3 fi?? 



9i?i „,(9i?: 



V ^ afc? ' ^fc? 



+ ^2rj,2 +^3;., 2 



dki 



(34) 
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Dividing the amplitude equations ( |T^-( [TTp through by i?i, -R2, -R3 respec- 
tively and then differentiating gives 



° - -^^^Wr'^'M^^'^'W'^'^'MJ'^^^'M 



dkl ^ ' dkl 



-27(^21^ + i?37T7| , (35) 



acosa (r, j,dRi dR^ dR2\ dR2 



-2,(i?3S + i?.||l, (36) 



dkl ^ ' dkl 



a cos a ( dR2 dRi dR^\ dRs 



dRi dR2 

dkl dkl 



-27(i?iini + ^2^ ' (37) 



which when multiplied by i?f, R2, Ri, respectively and added give 



I f)n dn d¥i 

= -{Rl + Rl) + acosa[^R2R, + Ri^R, + RiR2^ 



dkl ' ' ' 'dkl ■' ' ' ^dkl 



-27 U^-^^iRl + Rl) + R^^iRl + Rl) + +R^^(Rl + Rl 
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and hence 

dL 



dkl 
Similarly it is found that 

dL 

M 

dL 



(Ri + Ri). (39) 



(Rl + Rl), (40) 

2Rl (41) 



<9(ki.k2) 

and hence 

{Rl + Rl)Q,T + RlQ2T = V. (ki(i?2 + Rl) + k2i?2) , (42) 

{Rl + Rl)Q2T + RlQiT = V. (MRl + Rl) + k^Rl) , (43) 

which gives the Cross-Newell equations 

{RlRl + RlRl + RJRDe^T = RlV.{k^Rl-ksRl) + Rl\/.ik^Rl~k2Rl), 

(44) 
{RlRl + RlRl + RlRl)e2T = RlV.{k2Rl-ksRl) + RlV.{k2Rl-k^Rl) 

(45) 

upon rearrangement. 
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4 Cross-Newell equations for triangles 

A similar approach can be adopted to derive Cross-Newell equations for tri- 
angles starting from the governing equations 

— i = Zi{\ - ai\zi\'^ - a2{\zi\'^ + Izil"^ + \z^f) - a3{ziZ2Zs + zlz*z^)} 

+^:+i^l2{S - «3(ki|' + k2|' + ksH} + V=^^,. (46) 

These are the lowest order amplitude equations that permit triangles as a 



stationary solution Wi\ and once again the spatial derivatives have been 
chosen to ensure that the governing equations are isotropic. Fully nonlinear 
stationary triangles satisfy 

= X- kf - ailzil"^ - a2{\zi\'^ + \z2\'^ + Izsl"^) - as{ziZ2Zs + zlz2Z^){A7) 
= 5-asi\zi\^ + \z2\^ + \z3\^). (48) 

Writing Zi = RiC^^^ as before with a = Oi + 62 + 0^, gives 



. hM±M±M) + f , (49) 

Rl = S^M±M)+^, (50) 

6ai das 

Rl = M±|L:M + ^, (51) 

3ai 803 
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-^ - 6^^(^^-(*? + *^ + *^)-^'^)^ (^^' 



The Lyapunov functional is given by 

L = -X{\z,\'+\z,\'+\z3\') + ^{a^ + a2){\zr\' + \z2\' + \z3\') 

I n |2| |2 I I |2| |2 I I |2| |2\ /'cqN 

+a2(|zi| |2;2| +1^2! \zs\ +\zs\ \zi\ ) (53) 

I /I |2|| |2|| |2\/ I ***\ 

+03(^11 + F2I +\Z3\ ){ZiZ2Z3 + Z^Z^Z^) 

-6iz,Z2Zs + zlzlzl) + iV^ip + |V^2|' + |Vz3p, (54) 

which gives, upon substitution for Ri and a, 
1 



3ai 



{A;^ + A;^ + 4(ki.k2)' + klkl - 2A;2(ki.k2) - 2fc2(ki.k2)} 



25, ,2 ,2 , 1 X 35A 5(52(ai + 3a2) 

— (A;? + A;2 + ki.k2) + ^ 

as as bai 



+-{ki + ki + ki.k2) - — + ^"^":;^"^^ (55) 



As before the phase equations are given by 



2{Rf + RDOiT + 2^=92,- = V, ( 2kig^ + k^-g^^^^j , (86) 



2(Rl + Ri)e,r + 2flle.. = V. (2k.^ + k,^^) . (57) 



It is easily seen that 



^ = --^(2A:?-A:2^ + 2ki.k2) + -, (58) 

okf 3ai as 
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M - -3^(-*?^2^-5-*>-^^)^|. (««) 



^^ ^ ■(4k,.k2 + fc? + ki) + -. (60) 



t}(ki.k2) 3ai as 

After some rearrangements and substitutions the Cross-Newell equations for 
triangles are found to be 

i?2v.(ki(i?2 + 2(5/3a3) - MRl + 2^/303)) 

+i?^V.(ki(i?2 + 25/303) - MRl + 25/803)), (61) 

{R1R2 + -R2-^3 + -R3-Rl)02T = 
i?2V.(k2(i?2' + 25/303) - ki(i?2 + 25/303)) 

+i?2v.(k2(i?^ + 25/303) - ks{Rl + 25/303)). (62) 

5 Free hexagons and triangles 

The Cross-Newell equations take a different form when the total phase o is 
not constrained to remain fixed by the governing equations, for example in 
the case a = in §3. All three phases 6i are independent, with the following 
consequent modifications of the analysis 

Vx = kidg,+k2dg,+ksde,+eV^, (63) 
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dt = eieiTde, + e2Tdo, + e3Tde,), (64) 

Sz*Q = dejZ*Q59i + de2Z*Q592 + de^^ZgdOs, (65) 

which lead to the phase equations 

dG ' 



eiT|9,iu;o|2 = V. ( 2ki^^ , (66) 



dG 



&2T\de2Wo\' = V. ( 2k2T^ , (67) 



dG^ 



OsTldesWol' = V. ( 2k3Tn:2 • (68) 



In the hexagon case of §3, the Cross-Newell equations would be 

RlQiT = V.{Rlk^), (69) 

RlQ2T = V.{Rlk2), (70) 

i?2e3T = V.(i?2k3). (71) 

These are the equations that would have been found for three independent 
sets of rolls in the same system, as might have been expected, since nothing 
in the analysis constrains the hexagons or triangles to remain hexagonal 
or triangular. In particular, it is not to be expected that the condition 
ki + k2 + ks = will be maintained over long times. 

For hexagons which remain everywhere exactly hexagonal, such that 
|ki| = |k2| = \ks\ = k and Ri = R2 = R3 = R, the phase equations 



also take the form ( ^S|)-( |7T]), since the size and orientation of a hexagon can 
then be determined from a single wavevector, as in the roll case. However, 
this constrains the hexagons to behave as a rotating, shrinking or expanding 
lattice, which is clearly not a realistic model for most experiments. 

6 Flux-divergence form and the general non- 
variational case 

The Cross-Newell equations ( P) and ( ^ for hexagons and triangles in gra- 
dient systems are in flux- divergence form, which has consequences for defect 
formation, as in the case of rolls [§. Note that stationary solutions of equa- 
tions ( ^ and ( ^) take the form 

-■hi -"to) -■ (-) 

Following 0, it is interesting to set 

dG dG dG 



dkl dkl ^(ki.ka) 



1, (74) 



which implies that V.ki = V.k2 = hold. Since the wavevectors are gradi- 
ents of the phases, it is clear that Vxki = Vxk2 = also hold, and that 
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V^Gi = V^02 = 0. The solutions of these equations are the harmonic de- 
fects catalogued in p. Despite being energetically unreasonable, and hence 
looking somewhat unphysical, because they contain features at wavelengths 
which lie outside the stable region, they provide a good illustration of the 
topology of real defects. It is possible to construct a harmonic hepta-penta 
defect by positioning two harmonic dislocations |^ on top of each other, as 
shown in figure Q. The hepta-penta defect is the canonical defect of hexagons. 
In the general nonvariational case, the Cross-Newell equations can be 
written 

r(A;^,fc^,ki.k2)eiT = ai(fc2, fc^, ki.k2)V.ki + a2(A;^, fc^, ki.k2)ki.Vfci 

+a3{kl, kl ki.k2)ki.VA;2 + Mkl kl ki.k2)ki.V(ki.k2) 
+(3i{kl kl ki.k2)V.k2 + P2{kl kl ki.k2)k2.VA:i 

+(3s{kl kl ki.k2)k2.VA;2 -I- P^ikl kl ki.k2)k2.V(ki.k2), 

(75) 
T{klklk^.k2)Q2T = al{klklk^.k2)V.k2 + a2{klklk^.k2)k2.Vk2 

+a3{kl kl ki.k2)k2.VA;i + a^ikl kl ki.k2)k2.V(ki.k2) 
+(3,{kl kl ki.k2)V.ki + P2{kl kl ki.k2)ki.VA;2 

+(33{kl kl ki.k2)ki.Vfci + (3^{kl kl ki.k2)ki.V(ki.k2). 

(76) 
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In contrast to the variational case, these equations cannot in general be 
reduced to flux-divergence form, and hence it cannot be assumed that such 
general hexagonal patterns will have defects whose topology is given by that 
of harmonic defects. 

7 Conclusion 

This paper has derived the Cross- Newell equations for triangles and hexagons 
in a general real gradient system. The resulting equations can be put into 
flux-divergence form, indicating that the topology of defects of such a hexag- 
onal pattern can be described by that of harmonic defects [||. The general 
nonvariational case, however, is not flux- divergent. In both cases, the phase 
stability boundaries and modes of instability for general distorted hexagons 
and triangles can be recovered from the Cross-Newell equations. 

An explicit analytical form for the Cross-Newell equations is found for 
both hexagons and triangles in the case where the governing equations are 
generalisations of the corresponding complex amplitude equations. 

This work suggests avenues for further investigation. In particular, it 
would be interesting to analyse the Cross-Newell equations in a general non- 
variational system, and also to integrate the phase equations numerically in 
order to make a comparison with an integration of the full governing equa- 
tions, for example in order to compare the regions of stability of hexagons. 
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A further interesting possibility is to use the Cross-Newell equations to in- 



vestigate the simultaneous occurrence of up- and down-hexagons |T5[ . These 
avenues will form the basis of future work. 
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Figure 1: A harmonic hepta-penta defect. The figure shows contours of 
/(X) = cosGi + cos ©2 + cos 03, where 6i = -Rcos(q; + 7r/2) — (a + vr/2), 
62 = -Rcos(a + 7r/6) — (a + 7r/6), ©3 = —61 — ©2, and where (R, a) are 
polar coordinates for X. 
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